Excitation of magnetostatic spin waves in ferromagnetic films 



O 

o 

(N 
X 



Yu. Kuzovlev, N. Mezin and G. YaroslQ 
A.A.Galkin Physics and Technology Institute of NASU, ul. R.Luxemburg 72, 83114 Donetsk, Ukraine 

(Dated: February 2, 2008) 

The torque equation of nonlinear spin dynamics is considered in the magnetostatic approximation. 
In this framework, exact expressions for propagator of linear magnetostatic waves in ferromagnetic 
film between two antennas and corresponding mutual impedance of the antennas are derived, under 
conditions of uniform but arbitrarily oriented static magnetization and arbitrary anisotropy. The 
results imply also full description of spectrum of the waves. 
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I. INTRODUCTION 

The terms "magnetostatic oscillations" or "magneto- 
static waves" (MSW) 0, H H mean such relatively 
long spin waves (SW) whose properties are dominated 
by the long-range quasi-static dipole interaction between 
"spins" and eventually by geometry of ferromagnetic 
sample (its shape and dimensions). In opposite, proper- 
ties of relatively short SW (spin waves in narrow sense) 
[J 13 come rather from the short-range exchange in- 
teraction. In ferrite films, MSW can be easy excited, 
with the help of microwave-frequency magnetic fields in- 
duced by wire or strip- like antennas 0, 0, II] j up to the 
level of strongly nonlinear interactions of MSW between 
themselves and with SW. In addition, small velocity of 
MSW ensures compactness of nonlinear transformations 
of electromagnetic signals. When exploring these pos- 
sibilities, it is important to know characteristics of lin- 
ear MSW modes and corresponding linear impedances of 
their antennas. Concerning MSW in flat homogeneous 
films, there are widely used exact results by Damon and 
Eshbach Q ■ Other authors M, [ljj (see also paper by 
Kalinikos and coworkers in [f| and references therein) 
developed approximate approaches to films which are 
non-homogeneous because of specific surface effects. At 
the same time, to the best of our knowledge, the homo- 
geneous case remains incompletely investigated. In the 
present paper we find new solutions to this case believing 
that they will be useful for both practice and theoretical 
modeling of more complicated situations. 

We confine our consideration by "magnetostatic ap- 
proximation" which neglects SW (in narrow sense) at 
all, as if radius of the exchange interaction, r e , is equal 
to zero. This formal trick seems reasonable when cross- 
section dimensions of all the antennas or/and distances 
from them to film's surface are much greater than r e . In- 
deed, under this condition magnetic field created by an- 
tennas is so smoothly distributed in film's interior (with 
spatial scales much greater than r e ) that must only very 
poorly excite those SW whose length is comparable with 
r e or smaller. We do not know strict proof of this 



statement, but it rather convincingly follows from re- 
sults of the "magnetostatic approximation" itself (in par- 
ticular, results presented below) which demonstrate, at 
mentioned condition, insignificance of "infinitely short" 
MSW in the sense of both their private excitations and 
summary contribution to film's linear response to weak 
external field. Other deal is when in definite nonlinear 
processes even SW much shorter than r e can be generated 
from long MSW [|| (by the way, all the nonlinear models 
collected in ||| agree with the above statement). How- 
ever, here we are interested in linear (small-amplitude) 
response only, although examining the path to it from 
basic nonlinear equations. 

Besides, firstly, we will keep in mind that film under 
consideration is thick, in the sense that its thickness, D , 
is large as compared with r e . This condition will make 
our results be better applicable to real films with surface- 
induced non-uniformity. Secondly, if a ground (static) 
state of film's magnetization is divided into domains and 
hence strongly non-uniform then this state itself must be 
understood and described before its magnetostatic ex- 
citations. Of course, we avoid this extremely difficult 
problem, being satisfied by consideration of MSW on 
background of uniformly magnetized static state. For- 
tunately, for example, in yttrium-iron garnet films such 
state can be enforced by comparatively small bias field. 
At formulated assumptions, we will obtain the function 
("propagator of MSW") which describes mutual influ- 
ence of two antennas by means of MSW propagating in 
a film with arbitrarily directed ground magnetization in 
presence of arbitrary magnetic anisotropy. 



II. BASIC EQUATIONS 



The classical semi-phenomenological dynamics of mag- 
netization in solid ferromagnets and ferrites 0, 0, |n| 
is based on the Landau-Lifshitz-Gilbert equation (the 
torque equation) 



— = F x S 

dt 



7{F-S(S-F)} 



(1) 
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Here S is unit vector in the direction of local magneti- 
zation (\S\ = 1); F is thermodynamic force, or effective 
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internal magnetic field, expressed in units of the satura- 
tion magnetization, M s ; the time is expressed in units of 
To = (2%gM s )^ 1 (g w 2.8MHz/Oe is the giromagnetic 
ratio); x and • are symbols of vector and scalar prod- 
ucts, respectively; 7 is phenomenological friction (dissi- 
pation) parameter. For our purposes, the simplest model 
of dissipation is sufficient as those introduced by Eq^ 

The internal field, F, is composed at least from 
0, H, El (i) external bias field, H e , (ii) magnetic 
anisotropy field, H a , (iii) own magnetic field created by 
the magnetization (by "spins"), H s , and (iv) exchange 
field, H exc h ■ We do not add to this list magnetic field 
induced by electric currents in the ferromagnetic sample 
since assume that it is good isolator. Usually small di- 
mensions of real samples allow to neglect time delay of 
H s , therefore H s can be found from quasi-static ver- 
sion of the Maxwell equations. Particularly, in absence 
of conductors and other ferromagnets in vicinity of the 
sample, these equations yield 



GS = V V 



S{rl)drl 



r/\ 



(2) 



In most simple model of exchange interaction, H exc h = 
r 2 V 2 S 0,1a, 01- It should be noted that exchange in- 
teraction ensures local smoothness of the magnetization 
distribution, S , while the operator of dipole-dipole in- 
teraction, G , is bounded when acts on smooth distribu- 
tions: ||G|| < 4-7T . If A(S) denotes density of energy of 
the anisotropy then 



-A'(S) = -dA(S)/dS 



In principle, all the consideration in this section and 
next section can be easily generalized to non-uniform 
anisotropy if suppose that A(S) is some function of space 
coordinates. 

Let the subscript "0" be attribute of static magneti- 
zation state at constant bias field, H e = H a = const. 
According to Eq^in any such state the vectors So and 
i*b are parallel one to another, that is Fq = Wo So, where 
scalar field Wo (absolute value of static internal mag- 
netic field) is defined by the requirement |S| = 1. 

Let additional time-varying magnetic field is switched 
on: H e = Hq+ h, h = h(t) . Corresponding deviation 
of magnetization from its static value, s = S — Sq , can 
be represented as 



s = Sl + (S|| - 1)S , S± = nS _L S 



(3) 



n 



1 -So® So, S\\ = ±vi - \Si 



where (§) means tensor product, and II is operator (ma- 
trix) which performs projection of vectors onto the plane 
perpendicular to So- If none spins are overturned by the 
perturbation then sign of S11 is everywhere positive. In 
terms of s and Sj_ , Eq.l transforms into 



dS± 



S\\[F ± x So]+7(l-S ± ®S ± )F ± 



(4) 



where F± = — 5E±/5S±, and E± is the excess energy 
(energy of excitation) implied by the perturbation, 



E± = 



1 W s 2 + A(s) + ^-s-Gs-h-s + C(Vs) 



dr 
(5) 



2 " w 2 

Here function A is defined by 

A(s) = A(S Q + s)- A(S ) - s ■ A' (So) 



and function C represents exchange contribution to the 
excess energy. In the mentioned model, C(Vs) = 
^cEa^"^) 2 ' The functional derivative Fj_ in Eq0] 
should be evaluated with taking into account full depen- 
dence of s on Sj_ in accordance with Eq|3 

Importantly, the frictionless version of Eq0] (that is at 
7 = 0) follows from the variational principle 



So 



dSj_ 

dt 



X Si 



dr 



I + S11 



E±}dt = (6) 



Of course, in general So is a complicated function of 
spatial coordinates, hence all the related values (Wo, II, 
A(s), and so on) are space dependent. 



III. LINEAR WAVES 

If the functional © represents positively defined 
quadratic form then the static magnetization pattern So 
is stable with respect to any small perturbation, and 
therefore we can speak about linear eigenmodes of the 
excitation. In this case, let us introduce the spin preces- 
sion operator, i?, the anisotropy matrix, A, and besides 
the exchange operator, C , by the relations 

RV = S x V , CV ee — r 2 e V 2 V , 



A a f3 = d 2 A(S )/dS 0a dS 0l3 
For the linear regime Eq0| yields 
dS, 



dt 



= (R- 1 U)(WS ± -h) , 



(7) 



where we introduce the integral-differential operator 

Wee Wo + A + G + C 

Rejecting from Eq0 both dissipation and external 
pump, we obtain equations for SW and MSW eigenmodes 
and eigenfrequencies: 



Ve 



-iujV = RWV 



(8) 



Herewith it is sufficient to consider positive frequencies 
only. Let the eigenmodes be enumerated by an index 
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k . Since in a stable state the operator W is positively 
defined, we can write 



A 



o 



Wo 



to k V k = iW 1/2 RW 1/2 V k , V fc = W 1/2 V k 



(9) 



The operator on right-hand side in the left of these two 
equalities must be self-adjoint, hence, its eigenfunctions 
V k can be made mutually orthogonal. From here the 
orthogonality rule for the eigenmodes does follow: 



iJs -[V k x V; n ] dr = 5 mk 



(10) 



The same rule is specified by the variational principle Jfjjl . 

To get more general formulation of linear theory, we 
should return from the "ready" dipole interaction oper- 
ator G to Maxwell equations: 



dS± 
dt 



(R--yIl){(W + A + C)S±-h-h s } , (11) 



AirS i ) = , V x h s = 



(12) 



where we introduced new vector field, h$, which repre- 
sents time- varying part of magnetic field self-induced by 
magnetization (i.e. the same as the whole field induced 
by s). As before, it is assumed that the sample is non- 
conducting. 

Applying Fourier transform to Eci llll in the frequency 
representation (frequency domain) from l|ll|) and (|12ll we 

have 



Sx = x[h + h s ] , V • m [h + h s ] = , 



(13) 



x = + (R - 7n)(VFo + a + d)}-\R - 7 n) , (w) 

where ji = 1 + Att\ ■ Obviously, because of presence of 
the differential operator C in denominator of the polariz- 
ability matrix x m l|14|l . in fact x is an integral operator. 

At this point we go to the "magnetostatic approxima- 
tion" formulated and discussed in Sec.l. Concretely, we 
reject the exchange operator C from denominator of (11411 . 
Formally, this is equivalent to that the exchange radius 
r e turns into zero (of course, herewith we do not neglect 
exchange interaction since it remains responsible for the 
magnetization phenomenon itself). Strictly speaking, si- 
multaneously the static magnetization, i.e. the patterns 
So and Wo, also must be treated in this limit. But this 
does not matter in the case of uniform static magnetiza- 
tion which we will investigate below. 

After that, x turns into algebraic expression becom- 
ing literally matrix, and the problem reduces to purely 
differential equations for the field hs- Direct analytical 
calculation gives very simple expression for the polariz- 
ability: 



(Wp + Ax+ A 2 )U -A ± + iuR 
(W + A 1 )(W + A 2 )~^ 



(15) 



1+7 2 



Here A\ and A 2 are those two eigenvalues of matrix A± 
which correspond to the pair of its eigenvectors perpen- 
dicular to So and one to another: A±a\^ = A\p a\ t 2 ■ 
We enumerate them so that Sq ■ [cti x a 2 ] > . At prac- 
tically interesting ferrite samples 7 < 1CP 3 , therefore j 2 
plays no role. 



IV. PROPAGATOR OF MAGNETOSTATIC 
WAVES IN FILMS 



Let the time-varying field h is induced by some con- 
ductors which are placed outside the ferromagnetic sam- 
ple and carry a.c. currents /„ distributed with densities 
InJn (n — 1, 2, ...). Then we can write 

h = } h m I m , V • h n = , V x h n = — J n (16) 

Here h n is magnetic field created by unit-value current 
in n-th conductor. The same function determines volt- 
age (e.m.f.), e„ , induced in the n-th conductor by time- 
varying magnetization of the sample: 



(17) 



The fields h n as well as the self-induced field hs can be 
represented in the potential form. 

In the linear regime, when s — > Sj_, the response of 
the sample divides into sum of partial responses: 



h n 



-VU„, hs^J^h 



-VU S . 



where U sn is potential of the field induced by the sam- 
ple in response to influence by n-th conductor. After 
obtaining hs in company with S±_ from i|17|) we will 
determine mutual impedances of the conductors, Z nm , 
caused by their interaction through the ferromagnetic: 
e n = "^2 Z nrn I m (the hat marks convolution operators). 

Now concretize the sample as plate (film) whose in- 
plane dimensions much exceed its thickness, formally as 
infinite plate. At sufficiently large bias field, Ho, ferro- 
magnetic plate allows for stable state of uniform magne- 
tization. In real finite-size films, such the state is demag- 
netized at film's edges only, in strip-like regions whose 
width is few of D (D is thickness). This justifies the 
theory of plane MSW in infinite uniformly magnetized 
film. 

Let film is disposed in the region —D/2 < z < D/2. 
Naturally, make Fourier transform with respect to time 
and in-plane coordinates, x and y , marking transformed 
functions by tilde. Introduce designations k = {k x , k y } , 
V = {ik x ,ik y , V z } . In film's interior V 2 ?7„ = 0, there- 
fore the potentials of conductors have the form 



U n (k,z) = $ n (fe)exp{|fe|(a n z-^/2)} 



(18) 
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where a n = 1 (— 1), if n-th conductor is placed above (be- 
low) film, and form-factor describes distribution of 
n-th current. In combination with (|12|l and H17[l the lat- 
ter formula implies the relation between the impedances, 
from one hand, and values of the potentials taken at film's 
surfaces, from another hand: 



^nm — n / \k\U x . 
Z7T 



U!,k,cr n ^j U n ^-k,a n ^pj dk 

(19) 

Here dk = dk x dk y /(2ir) 2 (in contrast with U n , poten- 
tials U sn are frequency dependent). Thus, the potentials 
are taken at the surface most close to receiving antenna 
(n-th conductor). 

The Eq^H or equivalcntly, 



(V-$V)(U sn + U n )=0 



(20) 



should be solved under standard boundary conditions 
[3. Illj. To write the answer, introduce the unit-length 
vectors 

u={k x /\k\,k v /\k\,0}, *= {0,0,1}, 
and besides the matrix 



(21) 



As usually, the solution is composed by two exponents: 



M = 








v ■ jlv V 


fXZ 






Hzz 




~z ■ flv ~z 





U n = ^u n ±exp(q±z) , q± = X±\k\ 



A± = A ± A , A = 



A 




(22) 



(23) 



(24) 



It should be emphasized, however, that in general case (at 
arbitrary orientation of the vector So) the exponents q± 
are neither poorly imaginary nor poorly real but complex, 
that is MSW is not standing in Z-direction. 

If we took into account finite exchange radius r e and 
dealt with the operator- valued polarizability matrix (|14H , 
then in place of (|22|l we would get a sum of at least six 
terms, where transverse wave numbers of order of |fc| (as 
q± in 1)22(1 ) are more or less hybridized with real or imag- 
inary wave numbers of order of 7r/r e . From the point of 
view of our aims, such complication would be meaning- 
less accuracy. But it can be necessary when considering 
short SW or small-scale details of long MSW on back- 
ground of a nonuniform domain structure. Most natural 
approach to these tasks is direct analysis of the system 
of equations (fTTft and i(T2*|) . 

Consider the susceptibility matrix M. In standard 
spherical coordinates let 9 be the angle between Z-axis 



and Sq. In the AY-plane (film's plane), we introduce 
quantities and v± as cosine and sine, respectively, of 
the angle (counted clockwise) between projection of So 
onto this plane and the above defined unit vector u lying 
in it. Besides, consider the plane a\a^ perpendicular to 
So and, in this plane, define ip being the angle between 
plane ZSq and vector a± (definition of vectors a\ and 
a.2 was done at the end of Sec. 3). Further, introduce 
the quantities A± = (A 2 ± A\)/2. At last, let Vl be 
nominator of the polarizability matrix (|15f) . 

fi = {W + A x + A 2 )U - A ± + iuR 
In these designations 

fi zz = (W + A+ + A_ cos 20) sin 2 9 , (25) 

n„„ = ( W + A+ ) (v\ + i/jj cos 2 9) + (26) 

+A-{(vj cos 2 9 - v\ ) cos 2-0 - 2v\\y± sin 20 cos 9} , 

fi Z i/, vz = fix T wv± sin 9 , (27) 
fi x = sin 9 {A-v±_ sin 2ip - (28) 

- f|| (W + A + + A- cos 2-0) cos 9} 

These formulas make it evident that effects of anisotropy 
are determined by A- , while A + merely redefines the 
magnitude of static internal field, Wo- 

Then, for a given in-plane orientation of the wave, v , 
introduce characteristic frequencies by the following ex- 
pressions: 



loI = (W + A+) 2 -A 2 _ 



UJ 2 . = OJq + 47r0 2 



w i 2 =UJ l + 2n(£l zz + £l vv ) =p 



(29) 



(30) 



(31) 



T 2TtJ{n zz + n vv ) 2 - {2u x sin 9) 2 uj 2 , 



u>l = {lu 2 +uj 2 + (4.m>x sin 6») 2 }/2 



(32) 



The frequency u> u which is independent on the in-plane 
wave vector k is the uniform precession frequency. In 
terms of these frequencies, 



imu)v±_ sin 



A 



(33) 



(34) 
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Besides, below we will need in the determinant 

2iJ 2 — LU 2 — LU 2 , . 

A = det M = — 3 „ „ (35) 



It appears that A being quadratic function of the matrix 
elements of M and \ nevertheless always has simple pole 
only. 



J 



We omit trivial but tremendous evaluation of the sur- 
face potentials which appear in Eq^5] The result, for 
the surface closest to a given antenna, looks as 



U m (u, k, a n D/2) = $„(fc)P(w, k) 



(36) 



P(u,k) = 



1 



f-^vz 



1 + A + 2/x 2J .Acoth(A|fe|D) 
I 



Amuis i sin 9 



G(u,k) 



(37) 



where denominator in the latter expression is given by 

G(uj,k) = w§ - lu 2 + (uj 2 u - w 2 )Acoth(A|fc|D) (38) 

For brevity, we do mark dependencies of the factors A 
and A on oj and k as well as dependencies of ^1,2,3 on k 
(or, to be precise, on direction of the in-plane wave vector 
k). Combining these formulas and Ea ll9l for mutual 
impedance of two antennas located on one and the same 
side from film, we obtain: 



Z nm = 1 ^- I \k\$ n (-k)$ m (k)P(cj,k)dh 
2tt J 



(39) 



The latter formulas present main results of the paper 
and, as far as we know, can give useful addition to re- 
sults of Damon and Eshbach 0] and other authors (see 
Sec.l). Function P(u,k) is the required propagator of 
linear (weak) magnetostatic excitations from one antenna 
to another. At the same time, it contains complete in- 
formation about spectrum of MSW. The condition that 



J 



its denominator turns into zero, G(ui, k) — (in absence 
of dissipation, at 7 = 0) , yields a set of dispersion laws 
for all possible types of MSW. This will be the subject 
of separate work. 



V. MUTUAL IMPEDANCE OF WIRE 
ANTENNAS 

To be more concrete, consider relatively simple but 
practically interesting case of straight-line wire antennas 
which have round cross-sections and are parallel one to 
another and to film's surface. Besides, let they be ori- 
ented along Y-axis and located at A-positions x n , on 
the same side from the film and at distances p n from its 
closest surface. In this situation 

$„(fc) = (An 2 /ick x ) eyip(-\k x \p n - ik x x n )8(k y ) 

(here c is speed of light) , and Eg 1391 concretizes to 



Z nm [Ohm] 
v [cm] /[GHz] 



= 4ni / exp{-q(p n + p m )} 



(1 - A) cos(ga;) + (ji zu - p v 



r 



1 + A + 2p zz K coth(AgD) 



) sin(qx) dq 

q 



(40) 



Here on the left w is the film's width (formally infi- 
nite) measured in centimeters along antennas (i.e in in- 
direction), / is the frequency expressed in GHz, while 
on the right-hand side x = x n — x m , and the integral is 
taken over q = k x at k y — > 0. The latter means that 
matrix elements of the magnetic susceptibility matrix M 
and the functions A and A are calculated at ^ = {1,0,0}. 

Note that the dimensionless circular frequency to 
which enters all these functions is connected with fac- 
tual frequency / expressed in GHz by the relation (see 
Sec.2) 



<*> = ///(> ' h= (27TT )" 



gM s 



The impedance (|40|) possesses evident asymmetry with 
respect to sign of x, if (p ZL , — p uz ) 7^ 0. This is one 
more example of the non-reciprocity inherent to wave 
propagation in presence of static magnetization. 

One can see also that the impedance is a function of 
the dimensionless ratios D/p n and x/D only. Clearly, 
this is the consequence of specific scale invariance of the 
dipole interaction which is sensible to shape of the sample 
but not its size. 

In obvious way, one can generalize Ea l40l for many- 
element antennas each consisting of several (mutually 
parallel) wires with alternated signs of current in them. 
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It should be noted also that all the formulas (|25|) - (|35[1 will 
serve for analytical calculations (e.g. in next our publica- 
tion) , but if using computer it is sufficient to numerically 
calculate matrix (|21[) and then factors (|24[l and l|35l) only 
(by this reason we expressed the impedance in terms of 
these quantities). In general, of course, first of all one 

I 



Of course, here in the integrand v — {sign(q), 0, 0}, and 
Ao is defined in lEHll . 

For simple example, let us evaluate self-impedance, 
Z\\i of straight wire antenna (to be precise, the contri- 
bution to full self- impedance stipulated by the film), in 
the special case when bias magnetic field vector, Hq, lies 
in the film's plane. For concreteness, let it be oriented 
along y-axis. Besides, we assume that characteristic 
magnetic field of anisotropy is small in comparison with 
|i?o|+47rAf s , which allows to neglect effects of anisotropy. 
At last, if we direct the antenna in parallel to H then 
the impedance must be caused primarily by the so-called 
surface MSW discovered by Damon and Eshbach • Un- 
der above formulated conditions, in fc-plane these waves 
occupy the sector \k y /k x \ < y/4TrM s /\H \ ■ But, natu- 

I 



In this formula, i?n =RcZn is the film-induced con- 
tribution to resistance of the antenna [Z = R— 2nifL); 
magnetic field is dimensionless, i.e. expressed in units 
of M s ; the frequency belongs to the interval foUJ u < f 
< fou>t ; oj u is dimensionless frequency of uniform pre- 
cession and at the same time lower bound of spectrum 
of the Damon-Eshbach waves, while tot is upper bound 
of this spectrum (and spectrum of MSW at all) . Outside 
this frequency interval, there are no waves perpendicular 
to the antenna, and hence R\\ turns into zero (or, to be 
more precise, becomes comparatively small) . 

Let us notice that X — » 1 when / — > foOj t , therefore 
under condition p/2D < 1 the resistance (|43[) tends to 
infinity at upper edge of spectrum of the DE-waves. The 
matter is that here the group velocity of DE-waves, v g , 



must find the static magnetization vector, So, but this is 
also not a hard task for computer. 

For the case when two parallel wires are situated on 
the opposite parties from the film, evaluation of corre- 
sponding boundary potentials yields (in the same units): 



(41) 

I 

rally, sufficiently (infinitely) long antenna excites mainly 
the waves with \k y /k x \ — > which run perpendicularly 
to the field. Apparently, the latter case is the only case 
when the dispersion law of MSW (Damon-Eshbach waves 
with k _L .Ho, or DE-waves) can be written in the evident 
analytical form |2(: 



iu DE (k) = y/\H Q \(\H Q \ + 4tt) + 4^ 2 [1 - exp(-2D\k\)] 

(42) 

Here field and frequency are expressed in the dimension- 
less units introduced in Sec. 2. 

In this case there is a good analytical approximation 
for the integral (|47I)l which yields 



(43) 



I 

turns into zero, hence density of states (DE-wave modes) 
tends to infinity. 

Under the opposite condition, p/2D > 1 , this effect 
is canceled by sufficient weakness of excitation of short 
DE-waves. The presence in (|43H of the exponent which 
depends on geometric parameters of the system is even- 
tually consequence of scale invariance of the dipole inter- 
action. 

What is interesting, in the lower part of spectrum of 
DE-waves the resistance Rn is almost independent on 
the film's thickness, D, although seemingly the e.m.f. 
and thus the resistance must be proportional to amount 
of magnetic moments (spins) under excitation and thus 
to D . The matter is that the energy outflow from the 
antenna, p , is proportional to the group velocity, v g , of 



Z nm = ^ f°° e -\ q \ {Pn+ p m ) +iqx ( e -\ q \ D 2p tzz Aexp(X \q\D) \ dq 



wf ~ J_ x \ (1 + A) sinh(A|<j|D) + 2^ zz Acosh(A|gr|L>) J \q 



i?n[Ohm] 
w [cm] /[GHz] 



Airujt X ( 1 - X 



1-X 2 \1 + X 



p/2D 



In 



l + X 



1-X 



X 



(xL, 



/ 



Air 2 



w = 4- i ^ = vlHo|p/o|+47r) , uh 
Jo 



2tt 



7 



the DE- waves under excitation: p oc Dv g \S±\ 2 (where 
S± represents magnitude of spin precession). From the 
other hand, we can write p oc e\ jR\\ , while e\ oc D \S±\. 
Three above relations result in Rn oc D/v g . But, as it 
follows from (|42|l . group velocity of long DE- waves is pro- 
portional to the thickness, v g oc D. This is the reason 
for the indifference of R\i{u) — ► to u ) with respect to D. 

A simple analytical estimate for the inductance, L\\ 
= —ImZn/(2irf), can be deduced at p/D ~ 1 only, 
and then 

« -(2 + \H \M exp(-A)Ei(X) (44) 

w[cm| 

(Ei is the integral exponent function). Clearly, L\\ can 
be both positive and negative. 

VI. CONCLUSION 

In brief, we found (i) propagator l|37l) of magnetostatic 
waves (MSW) running in infinite ferromagnetic film from 
one antenna to another and (ii) linear (small-amplitude) 
mutual impedance of the antennas, under arbitrary ori- 
entation of uniform static magnetization of the film and 
arbitrary magnetic anisotropy. 

Additionally, the equation G(u>, k) — with G being 
the denominator (|38|l in (|37() determines dispersion laws 



for various linear eigenmodes of MSW thus allowing gen- 
eralizations of classical results obtained in 0. 

To conclude, let us touch the applicability of formulas, 
obtained for continuous spectrum of MSW in infinite film, 
to real finite-size films where MSW spectrum is discrete. 
When in-plane dimensions of a film decrease then charac- 
teristic frequency separation between neighbouring MSW 
modes increases, but, at the same time, selection of the 
modes by any simple (for instance, straight-line wire) an- 
tenna becomes more and more worsened. As the result, 
the formulas derived for infinite film can give good es- 
timate for impedances of antennas interacting with real 
films. 

This expectation was confirmed by comparison be- 
tween the analytical estimates and measurements of 
impedances induced by millimeter-size ferrite films as 
well as by results of their numeric simulations. More- 
over, numeric simulations of the torque equation 
with dipole-dipole interactions between spins show that 
spatial-temporal patterns of spin precession even in 
rather small films (with length to thickness ratio ~ 30 
-j-100) and even at essentially non-linear regimes possess 
clear imprints of qualitative and quantitative character- 
istics inherent to linear MSW modes in infinite system. 
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